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Abstract 



A far-wing line shape theory based on the binary collision and quasistatic approxima- 
tions that is applicable for both the low- and high-frequency wings of the vibration-rota- 
tional bands has been developed. This theory has been applied in order to calculate the 
frequency and temperature dependence of the continuous absorption coefficient for fre- 
quencies up to 10,000 cm -1 for pure H 2 0 and for H 2 0-N 2 mixtures. The calculations 
were made assuming an interaction potential consisting of an isotropic Lennard-Jones part 
with two parameters that are consistent with values obtained from other data, and the 
leading long-range anisotropic part, together with the measured line strengths and transi- 
tion frequencies. The results, obtained without the introduction of adjustable parameters, 
compare well with the existing laboratory data, both in magnitude and in temperature 
dependence. This leads us to the conclusion that the water continuum can be explained in 
terms of far-wing absorption. Current work in progress to extend the theory and to validate 
the theoretically calculated continuum will be discussed briefly. 


i 



1. Introduction 


The absorption of radiation by water vapor in the window regions of the Earth’s 
atmosphere is a problem of great practical as well as theoretical interest (Deepak 
et al., 1980). In the past few years there have been a number of general reviews 
(Hinderling et al., 1987; Varanasi, 1988a; Grant, 1990; Thomas, 1990) concen- 
trating primarily on the experimental characterization of this absorption; for his- 
torical reasons this absorption is called the “water continuum” and the interested 
reader is referred to these reviews for original literature citations. As a result of 
this extensive work, there is nearly unanimous agreement on the density depen- 
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dence (quadratic), general agreement on the temperature dependence (strong, 
negative, at least in the 8-13 fim window for atmospheric temperatures), but 
considerable disagreement as to the magnitude and physical mechanism respon- 
sible for the absorption. Water dimers (Penner and Varanasi, 1967), polymers 
(Carlon, 1981 ), collision-induced absorption (Bimbaum, 1985), and the super- 
position of the far wings of collisionally broadened allowed dipole lines (Elsasser, 
1 938 ) have all been suggested as possible sources of the water continuum. While 
the first two mechanisms are undoubtedly present and contribute to the absorp- 
tion in the atmosphere (in varying amounts in the different spectral regions), 
there is increasing evidence to believe that they are not the primary mechanism. 
Although collision-induced absorption arising from binary collisions between 
water molecules would have the correct quadratic density dependence and strong, 
negative temperature variation, its magnitude is too small to explain the ob- 
served absorption. The dimer model, on the other hand, has a temperature de- 
pendence similar to that observed near room temperature in the 8-13 /mi win- 
dow (Roberts et al., 1976; Varanasi, 1988a; Thomas, 1990). However, this 
predicted temperature dependence is not in agreement with that of the observed 
continuum in other spectral regions (Burch, 1982; Thomas, 1990). Furthermore, 
the density dependence in the wings of the bands (presumably responsible for the 
continuum absorption ) would vary as the cube of the density in disagreement 
with the experimental results; also no spectral features unambiguously assigned 
to dimers have been identified in the laboratory data or in atmospheric 
measurements. 

Within the past decade, substantial progress for calculating the absorption far 
from the line centers of resonant transitions has been made. Starting from the 
basic formalism of Davies et al. (1982), Rosenkranz (1985, 1987) has devel- 
oped a quasistatic theory that is applicable in the high-frequency wing of the pure 
rotational band (the 8-13 fim window). This theory successfully predicted both 
the magnitude and the strong, negative temperature dependence. This theory was 
extended by Ma and Tipping (1990a), but, because of the nature of the approx- 
imations made, was limited to the high-frequency wings of the pure rotational 
and/or the vibrational bands. Utilizing the same basic formalism, but proceeding 
in a different way mathematically, these authors also developed a theory (Ma and 
Tipping, 1990b) that was applicable for the low-frequency wing of the pure ro- 
tational band, i.e., for the millimeter spectral region. Again, good agreement in 
both magnitude and temperature dependence between theory and experiment 
(Liebe, 1989) was obtained. 

More recently, a new far-wing line shape theory based on the binary collision 
and quasistatic approximations was formulated (Ma and Tipping, 1991) and 
applied to calculate the absorption coefficient, a(co), for frequencies to between 
300 and 1 100 cm - 1 for several temperatures. The results were in excellent agree- 
ment with the extensive experimental data of Burch and co-workers ( 1 984, 1 985 ) . 
In particular, by treating the resonant and anti-resonant terms separately, better 
agreement with the high-temperature laboratory data was obtained as compared 
to the theory of Rosenkranz (1985, 1987). This new theory has been generalized 
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to vibration-rotational bands (Ma and Tipping, 1992a) and to foreign broaden- 
ing (Ma and Tipping, 1992b), and thus constitutes a unified theoretical frame- 
work within which one can calculate the magnitude and temperature dependence 
of the self and foreign broadened far-wing absorption. 

Prior to the work of Rosenkranz (1985, 1987) and the development of the 
present theory, several previous attempts to model the continuum absorption in 
H 2 0 and H 2 0-N 2 mixtures have been made. Roberts et al. (1976) derived an 
analytic expression for the continuum in the 8-13 //m region by fitting the data 
of Burch et al. (1971, 1974, 1975). This formulation is still widely used because 
of its simplicity, but it does not reproduce the improved experimental laboratory 
data of Burch and co-workers (1982, 1984, 1985), nor does it model correctly 
the temperature dependence observed by Montgomery ( 1978) and by Loper et 
al. (1983). Thomas and Nordstrom ( 1980, 1982a, b, 1985) sought to model the 
continuum absorption in terms of a modified line shape. They assumed a Loren- 
tzian line shape near the line center matched to a statistical line shape in the far 
wings using a large number of adjustable parameters fit to various spectral re- 
gions at different temperatures. When extrapolated to other regions or to differ- 
ent temperatures, their results are not in agreement with experimental data. A 
more sophisticated line shape model was developed by Clough and co-workers 
( 1 980, 1983) who attempted to include as much physics as possible ( fluctuation- 
dissipation theorem, Nyquist theorem, etc. ) in their model. More recently, based 
on the early theoretical work of Davies et al. ( 1 982 ), Clough et al. ( 1 989 ) intro- 
duced an empirical line shape correction factor ( x function) that can be used to 
calculate the continuum contribution within the framework of FASCODE line- 
by-line calculations. This model, referred to as the CKD model, is based primar- 
ily on the data of Burch and co-workers, but over the years has been refined as a 
result of atmospheric comparisons. While it is certainly true that this approach is 
at present the best way to obtain the most accurate results over a wide range of 
frequencies without the introduction of a large number of ad hoc parameters, 
there are a number of drawbacks associated with this model. By introducing a 
single x function that is symmetric in frequency displacement, the resulting ab- 
sorption coefficient does not satisfy the principle of detailed balance (Ma and 
Tipping, 1994); this will hot affect the results in the near-wing region, but can 
affect them in the far-wing region. The temperature dependence of the empirical 
X function was obtained by fitting laboratory data at two temperatures (296 and 
338 K); these results can be interpolated successfully in temperature over this 
narrow range, but the extrapolation to the much lower temperatures in the at- 
mosphere can introduce additional uncertainties. In fact, there is some evidence 
(Kilsby et al., 1992) that the temperature dependence in the CKD model is re- 
sponsible for the underestimate of the absorption in the 10-12 /im window by as 
much as 30% in a tropical atmosphere. Extrapolation of the model to high-tem- 
perature data can also lead to difficulties (Hartmann et al., 1993). 

In the present paper, we do not make any attempt to compare the continuum 
models or the results obtained by the various authors mentioned above with our 
current results; neither do we compare our results with predictions made based 
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on alternative mechanisms (e.g., the dimer hypothesis). Rather, our purpose is 
to present the most recent theoretical results and to compare these with extensive 
experimental data in order to investigate whether far-wing absorption can ac- 
count for the major features of the observed continuum absorption. Accordingly, 
the outline of the present paper is the following. In the next section, we first re- 
view the mathematical formalism and discuss the nature of the approximations 
made. In Section 3, we then discuss the computational details and present theo- 
retical results for the positive and negative resonance theoretical line shapes and 
the resulting absorption coefficient a{u>) up to 10,000 cm" 1 for self and N 2 
broadening. In Section 4, we compare our present theoretical results with selected 
experimental data. In this regard, it is important to keep in mind that we have 
not made any attempt to improve the agreement of our results with specific ex- 
perimental data by adjusting the values of the input data. We also discuss briefly 
some preliminary validations of our theoretical results using atmospheric data, 
and conclude by mentioning some recent progress made to refine the present 
theory. 


2. The general formalism 

2. 1. The absorption coefficient and spectral density 

As is well known, if we divide a gaseous sample into an absorber molecule a 
and treat the remaining molecules as a bath, then the absorption of radiation of 
frequency cofcm" 1 ) per unit volume of a sample in thermal equilibrium at tem- 
perature Tis characterized by the absorption coefficient a {(d): 

4tt 2 

(x(a>) =-j-n a oj tanh{hw/2kT) [F(a>) +F{ — co) ], ( 1 ) 

5nc 

where c is the speed of light, « a is the number density of absorbers and F(a>), the 
spectral density, is given by 

oo 

F{aj)=-ReTr f [p a -e'^- L u pp a ]dt. (2) 

7t J 

0 

In this expression p a is the dipole moment operator of the absorber molecule 
and L is the total Hermitian Liouville operator which can be expressed as the 
sum of its components 

F=L 0 + Z., = L w + L w +L l . 

We assume 


Then Eq. (2 ) can be approximated by 
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OO 

F(co)=-ReTr f [p a -e i(w ~ Lo), e~' L ' , pp a ]dt. 
n J 


(4) 


By making a Laplace transformation of the operator e ,Ul pp : 


e~' Lu pn. 


-- f 

a 2 ni J 


- ico' t_ 


l 


oj' — L, 

— OO 

Eq. (4) can be rewritten as 


PPadof, 


(5) 


oo 

F(«)=^Im~ f d co'Tr[p a - ■ - — \ - — X — y p p a ] ■ 
7T Z7T7 J OJ — CO — L 0 CO — L\ 


( 6 ) 


We assume that the total density operator can be approximated by the product 
of three parts 

p = p^)p(b) pW ^ ( 7 ) 

where /> <a) and p (b) are the parts related to the unperturbed absorber and bath 
molecules, respectively, and p (ini) is the part related to the interaction between 
them. As a result we obtain an important expression for the spectral density in 
which formally all the variables except the internal variables of the absorber have 
been isolated 


oo 

F(m)=-Im- 1 -. f d oj' T r[p a -X 

n 2 m J oj-oj - L a 


<-r L r P 0nt) > b p (a) k)a. (8) 

The Liouville operator < ( 1 )/(co' — L, )/7 (int) > b contains all the information 
about the bath molecules and the interaction dynamics between the absorber 
molecule and the bath. Formally, this expression for the spectral density contains 
an integral. However, the integral does not complicate the calculations since it is 
easily carried out using the Cauchy integral formula. 

Based on the binary collision approximation which is valid in the low-density 
limit, we can focus on the absorber molecule and one bath molecule. Then, in Eq. 
(8) and subsequent expressions, the Liouville operator L { and the summation 
index b are only related to these two molecules. 

Since the isotropic interaction does not depend on the rotational quantum 
numbers, the Liouville operator associated with it will cancel out when it acts on 
the line space. Therefore, F is (r) need not be included in the interaction Liouville 
operator L,. The Liouville operator corresponding to the anisotropic interaction 
between the absorber and the bath molecule is 
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L i = [V, n M)I*-lV: nis (?)]/fi ( 9 ) 

= (G ams I*-IG* anis )/r\ 

where / (=/*) is the identity operator and the exponent s is the power of the 
leading term in the expansion of the anisotropic interaction. For example, s= 3 
for the dipole-dipole interaction (H 2 0-H 2 0), 5=4 for the dipole-quadrupole 
interaction (H 2 0-N 2 ), etc. In the following, We denote by |«> and G a , the ei- 
genvectors and eigenvalues, respectively, of the matrix operator G anis . 

The resolvent operator < ( 1 )/(to' — L, )/? (in,) > b in Eq. (8) is an operator act- 
ing on the line space of the absorber molecule a. Therefore, in general, it can be 
expressed as 


co—L x 


-p (inl) >b= Z 

ii jinj\ 


< — l ~rP (in,) > 6 1 i\r? > I hA > < i\r? I , (10) 

where |z, jj > are the basis in the line space of the absorber molecule, and i x 
stands for all quantum numbers specifying its states. By explicitly writing the 
summation over the internal variables of the bath molecule, its components can 
be written as 


< UA I < 7 ~rP (in,) > * I i\j? > = Z Z <iih\axfi\jih 

tO — L i i 2 j 2 a (I 

> t\ 71^“ m) l < « I * I h > <j\ h I /?> gj2 e- E * /kT /Q b , (11) 

\ja— {j a p/r J r 

where G alj =G a — Gp, i 2 and j 2 stand for all the quantum numbers of the bath 
molecule specifying its states, g j2 are the degeneracy factors due to nuclear spin, 
E j2 are the energies associated with the state \j 2 > , and Q b is the partition function 
for the bath molecules. The remaining summation over the external variables of 
the molecules is indicated by the trace with an index r. We note that the summa- 
tion over the internal bath variables in the internal bath line space in Eq. (11) 
includes only the states | i 2 i\ > because of the assumption that the bath molecule 
is not involved in the transition. 


2.2. The quasistatic approximation 

We assume that the interaction between the absorber molecule and a bath mol- 
ecule, V(r), is 

V(r) = V anis (r) + V is (r), (12) 

where F ams (F) is the leading term of the expansion of the anisotropic potential 
and V is (r) is a simple Lennard-Jones type model 


R.H. Tipping, Q. Ma /Atmospheric Research 36 (1995) 69-94 


75 


F is (r) = Ca- 6 [(a/r)'+( C 7/r) 6 ], (13) 

where C, a and t are parameters which are consistent with values obtained from 
other data (Rosenkranz, 1985, 1987). 

Based on the statistical theory in which the duration of a collision is assumed 
to be infinite and, therefore, which is valid in the far-wing limit, the relative mo- 
tion of the two molecules is not treated quantum mechanically, rather one uses 
the quasistatic approximation. In this approximation, the trace with the index r 
in Eq. (11) is replaced by an integration over r, the separation between two in- 
teracting molecules, and the part of the density matrix related to the interaction, 
Pa M) , is a statistical weight for a pair of molecules (the absorber and the bath 
molecule b) at the separation r and with their orientations specified by or, this is 
given by 

p( l inl '>=e- Vaa(r)/kT /U, 

where 

y aa (r) = <a | Ta nis (r)|a> + F is (r) 

= G a /r*+V ls (r). 

The normalization factor U is given by 

U=Tr[e~ v “°< r)/kT ] ra 

''max 

= J e - ViF)/kT dr. 

0 


(14) 

(15) 

(16) 


In the last step in deriving Eq. ( 1 6 ) , we have used the fact that the trace over 
rotational states related to the relative motion of the two molecules corresponds 
to a classical average over the angular variables of the relative motion. This 
expression can be rewritten in the form 


U=— J (1— e V(r>/kT )dr+ J dr 
0 0 
'max 

= - J (l-e~ v ^ /kT )d r+V, 

0 


(17) 


where V is the volume which can be replaced by 1 /n b , i.e., by the inverse of the 
number density of bath molecules. Thus for moderate densities and tempera- 
tures, U^l/n b and we obtain 

p< in "=n b e- v ^ /kT . (18) 

Therefore, within the quasistatic approximation, we have 
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CO 

? J- r; ^ int) l =4nn h \r 2 dr 

L o)-G a0 /r s X J « 


1 


,-Vaa(r)/kT 


o)-G a Jr 


( 19 ) 


With Eq. (19), the components of the resolvent operator <(1)/ 
(co— L, )p (,n °> b can be rewritten as 


<ij\\< — [ -y~P (int ] > b I i \j'\ > = 47t« 6 f r 2 dr £ £ < i , i 2 \ a > 

CO Lr\ J i 2 j2 dp 

<p\j\h > 1 j ,-v« a (r)/kT x < a | i\ j 2 > <j\ j 2 \P> gj2 e ~ E]2lkT l Q b ■ ( 20 ) 

It is easy to obtain their imaginary parts with the well-known formal identity 
lim ^-r^=P— ^ — inS(a- L,), (21) 

e^oCD— L,+ie CD— L, 

where P indicates the Cauchy principal value of the integral and <5 is a Dirac delta 
function having the Liouville operator cd—L, as its argument; explicitly 

Im<i y j\ | < 1 p u ’" ) >b\i\j?> (22) 

(O — L j 

GO 

= -4 n 2 n b fr 2 dr£ '£<i l i 2 \a> <0\ j l i 2 >S(co-G a0 /r s )e- l ' aa(r)/kT 

“ nh 

X < a | i\ h > <j\ j 2 1 P> g j2 e- Ej2/kT /Q b 

oo 

= -4n 2 n b \r 2 drY j £ < h h I a > < P\Ji h 

• nh afi 

> ^^S[r-{G afi /coy' s ) e - y ^ kT 
X <a\i\j 2 > <j\j 2 \P>gj 2 e~ E,2lkT /Q b 

4 7l 2 1 

= ——«*-£ X<i'|/ 2 |a> <P\j\h 

S COi 2 j 2 exp 

> ( \ G afS \/(o) 3/s e~ Vaa{( ' Gali/w) ' ls]/kT 
X &( G a p/ w)<a\i\j 2 > <j\ j 2 \P> gji e ~ K,2/kT /Q b , 
where &(x) is the unit step function. 
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2.3. The individual line shape functions 


We introduce the individual line coupling functions, Xtvivnn ( w )> defined by 
471^ 

Xuvnn ( ) — r -«* I Z #( )cu 3)/i | | 3/5 < / , i' 2 1 a > 

5 1272 «/? 

</J|7.«2> xe- l '°°" G °"“ )l/s '' kT <a\i\j 2 > <j\h\p>g j 2 e- E »' kT /Q b . ( 23 ) 

Then the imaginary components of the resolvent operator < ( 1 ) / 
( a )' — L, )/j (,n,) > b can be simply expressed as 

Im<iJ} | < 1 p m) > b \i\j?> = XiumA(°>)- ( 24 ) 

CO — L| CO 

With the Cauchy integral formula 

(25) 

it is easy to carry out the integration in the spectral density and obtain 


oo 

F{co)=-Im^—. [ do)’Tr[p a \ — — X 

n 2 ni J o)—(o —L„ 


■p 0M) > b p (a) fl a ]a 


(o’ —L i 


oo 

= d(0 '^ £ <>' l/Cl*l> y - - d m < ij T 1 

71 271/ J || yj CO CO CO I Jjj 


/> (in,) >»|/'i77>xa 1 e- |S,/ * 7 '<r, |£|/, 


o)—L l 


OO 

<2 ' lAi<1 

— OO 

X&, e~ m,kT <i\ \p a \j\ >/Qa=\x Z < 7 i l/C I b > 777 !. y fe im 


(w-m, Vl )g„ , c E *' kT < i\ | Ha | j\ > /Q a . (26) 

We note that the arguments of each line coupling function in Eq. (26 ) take the 
desirable form after the Cauchy integration over co' is carried out. 

More conveniently, the spectral density can also be expressed in terms of the 
individual line shape functions as 
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?((*)=]_ Z \ 77 ~ r r#iAu-<0i j) 

7l(Oij >0 tQij ) 

+e ^ ,kT T ^^ i{t0+(0 A Pi 1 ^ 1 2 ' (27) 

Here the individual line shape functions are 
47 r 2 

XiyjA tO)=—n b Y, E Z Z Z( <a l i l i 2> <*1 \Pm\j\ > 

S m mi x m j\ &P f 2 


<;./2U?>)*6»(^/ tu )w <i - 3)/i |G^| (28) 

Xe- VaaHGa ^ /o>) ' ls]/kT ca\p (b) p (a) fi m \p> / p u \p iv , | 2 , 

where m h and m jl are the magnetic quantum numbers related to the state i and j, 
respectively; for simplicity we have omitted the subscript 1 which indicates the 
absorber molecule. If we define the symmetric individual line shape functions 


Xim^)^ ,2kT XnAoj), 

then the absorption coefficient «(&)) can be written in the familiar form 

, , ^ ^ (osh(fiu>/2kT) if 1 

a(co) = n a Z Su T7I — ~ V 7 . ^ (m-m,,) 

oj~to (DjjShifuDij, !2kT) n ((co-ci),j) J 1 


(29) 


[ 

((O + CO'j) 


2 xAco+Wij) >, 


(30) 


where S v are the line strengths. 

We note that the individual line shape functions are not even functions of their 
argument 

tu (o})=Xji(-0)) but Xij(o9)^Xiji-co). (31) 

implying that the line shapes are asymmetric about their centers. 


2.4. The band independence approximation 


In the present paper, we are interested in the case in which the absorber mole- 
cule has a large permanent dipole moment and the leading term in the anisotropic 
interaction is due to the interaction between this dipole and the leading multipole 
moment of the bath molecule. In light of the fact that the matrix elements of the 
dipole moment of the absorber molecule between different vibrational states are 
at least an order of magnitude smaller that those between the same vibrational 
states, and that the matrix elements of the dipole moment between the same vi- 
brational states are only slightly dependent on its vibrational quantum number, 
we can assume that the matrix of G takes a block diagonal form. Using this as- 
sumption, the lines belonging to the different vibrational bands are not coupled 
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FREQUENCY a (car') 

Fig. 1 . The positive and negative frequency resonance average line shape functions for H 2 0-H 2 0 self- 
broadening (in units of cm - '/atm ):*+(«) (dashed curve) and ;£_(«) (dotted curve), as a function 
of frequency w ( in units of cm “ 1 ) calculated for T = 296 K. The solid curve is the intensity-weighted 
average line shape function. 

to each other, and the lines within each band are coupled in the same way. More 
precisely, the latter means that the individual line shape functions are common 
for all the rotational lines in a given band. If we denote those obtained from the 
pure rotational band by ^ y (oj) where the indices i and j indicate the rotational 
quantum numbers only and indicate the vibrational quantum numbers by v, and 
Vj, then the general expression of the individual line shape functions can be writ- 
ten as 

Xvuvjiioi) =Xvi,iVj, j {(o)=X ij {(o) (32) 

Based on the argument above, the total spectral density can be expressed as the 
sum of its components corresponding to each individual band 

F(co)='£F ViVj ((o). (33) 

Vi Vj 

Furthermore, the components of the spectral density, F ViVj (co), can be ex- 
pressed as 
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Fig. 2. The positive and negative frequency resonance average line shape functions for H 2 0-N 2 for- 
eign-broadening (in units of cm _l /atm): x + (w) (dashed curve) and x_(w) (dotted curve), as a 
function of frequency t o ( in units of cm _ 1 ) calculated for T— 296 K. The solid curve is the intensity- 
weighted average line shape function. 


FviVj{(0)— X \(. \23!iy( w Wvii.vjj) 
ft to Vii, Vjj > 0 UJ Vji, Vjj ) 


-J- e f>ajy,i,Vjj/kT_ 


(OJ+CO VllVjJ ) 


2 Xji(W+(0 Vii yj) rp Vii I Pv,, 


KVj) I 


( 34 ) 


Finally, the total absorption coefficient can also be expressed as the superposi- 
tion of the components corresponding to each individual band; i.e. 


a(co)= £ avivj(co), 

ViVj 


( 35 ) 


where 
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FREQUENCY w (cm' 1 ) 

Fig. 3. The H 2 0-H 2 0 absorption coefficient a(w) (in units of cm 2 molecule -1 atm -1 ) as a function 
of frequency a) ( in units of cm - 1 ) calculated for 7'= 296 K. The results obtained from the two aver- 
aged line shape functions are the solid curve, and the results from the intensity-weighted average line 
shape function are the dotted curve. 



FREQUENCY « (cm- 1 ) 

Fig. 4. The H 2 0-N 2 absorption coefficient a (w) (in units of cm 2 molecule -1 atm -1 ) as a function 
of frequency cu (in units of cm - 1 ) calculated for T = 296 K. The results obtained from the two aver- 
aged line shape functions are the solid curve, and the results from the intensity-weighted average line 
shape function are the dotted curve. 
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Fig. 5. The Hj 0-H 2 0 absorption coefficient a(w) (in units of cm 2 molecule -1 atm -1 ) as a function 
of frequency co (in units of cm -1 ) calculated for T= 296 K; the experimental values of Burch et al. 
( 1 984, 1985) are denoted by +, while d corresponds to the theoretical values calculated with two line 
shape functions, and O to those calculated with one line shape function. 


. y 9 . cosh(hco/2kT) 1 f 1 

oty iV j{o) n a ^jL >0 Ki,> Vv (o Vii yj sh ( ha)y iit VjJ / 2kT) n\(co-o) ViitVjj ) 2 

*</(«- Q) Vli , Vji ) + 7 —— nXjAaJ+oj^^) [. (36) 

In general, the line shape functions decrease very rapidly when their arguments 
increase; therefore, at the frequencies of interest the main contributions to the 
total absorption coefficient come from the nearby strong bands only. This is es- 
pecially true for the cases in which the leading anisotropic interaction between 
two interacting molecules is weak. 

In practical calculations, we can also introduce the positive and negative fre- 
quency resonance-average line shape functions (Ma and Tipping, 1991, 1992a, 
b), x + (co) and X- (m), which are common for every band and which simplify 
the calculation of the absorption coefficient. We note that the first terms in the 
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FREQUENCY u (cm' 1 ) 

Fig. 6. The H 2 0-N 2 absorption coefficient a(o>) (in units of cm 2 molecule -1 atm -1 ) as a function 
of frequency w (in units of cm -1 ) calculated for T= 296 K; the experimental values of Burch et al. 
( 1 984, 1985) are denoted by + , while A corresponds to the theoretical values calculated with two line 
shape functions, and O to those calculated with one line shape function. 


curly bracket in Eq. ( 36 ) are related to the positive frequency resonances and the 
second terms to the negative frequency resonance. We define the intensity- 
weighted positive frequency resonance-average line shape function by 

*+(& ))= I Xu(ai)p, I Hu 1 7 £ Pi\<i\M m \j>\ 2 (37) 

coij > 0 ij 

Et > Ej 


=— I I I E(<alM2> <h\p m \ji > 

S iiji m gl{ 1 i2 

Eii >Ej t 

<hh \P> )*0(G a »w (s - 3)/s |G^ I 

■^g — fico[(Ga + Gp)HGa — Gp)]/2kT— Visl(Gafi/co) l ti]/kT 

X<a\p w p (a) (i m \p>/ £ Pt\<i\Pm\j>\ 2 - 

u 

Ei > Ej 
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Fig. 7. Comparison between the H 2 O-H 2 0 theoretical values and the laboratory experimental data of 
Burch et al. (1984, 1985) for three temperatures in the spectral range 2400 cm -1 <w<2700 cm -1 . 
The solid curves are the theoretical results for temperatures T= 296 K, 328 K and 428 K, from the 
top to the bottom, respectively. The corresponding experimental data are denoted by 1, x, and J, 
respectively. 

Similarly, we define the intensity-weighted negative frequency resonance-av- 
erage line shape function by 

X-((o)= X Xji(o>)Pj\l3j,\ 2 / X Pj\ <j\Mm\i> I 2 (38) 

coij > 0 ij 

Ei > Ej 

drr ^ 

= n b X X I X( < «l'l I 2 >< i. \Pm\jl> 

S i\j\ m ap i*2 
Ei\ cEj t 

<hy 2 \P>)*0(G ap /(0)^ s -^\G afi \ 

y^ e -flw[(Ga+Gf)/(Ga-Gf)]/2kT-yisUG a ll/w)U-i]/kT 

X<a\p (b) p (a) n m \P>/ X Pi\ <ilA*«lj> I 2 - 

ij 

Ei<Ej 
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FREQUENCY ufm' 1 ) 

Fig. 8. Comparison between the H 2 0-H 2 0 theoretical values and the laboratory experimental data of 
Burch et al. (1984, 1 98 5 ) for T= 296 K. in the spectral range 3000 cm - 1 < co < 4300 cm ~ 1 . The theo- 
retical values are denoted by A and the experimental values by +, respectively. 

We need only to evaluate x+(co) and %_ (co) for a» 0 since the values of these 
functions for w<0 are approximate equal or equal to their partner’s value for 
Gt>0; i.e., /+(-") =*-(")• If we do not distinguish between these two aver- 
ages, we define the intensity-weighted line shape function /(co) 

_ 47 z 2 

%((0)= 'Z.Xij ( CO )Pi I flij 1 7 Z A I < « I L/> 1 2 = —n b Z < a I Pm I P 

ij ij A cifi 

>*&(G /i /(O)C0 (S ~ 3)/s \G ^ | ^CSy^g-fitvUGa + Gfi)/ (Ga — Gp) ]/2kT— Vj S [(G a p/ oi )> / 3 ] / kT ^ 

<a\p w p {a) ii m \p> /'£p,\<i\Pm \j> I 2 - 

(39) 


With the two average line shape functions, the components of the absorption 
coefficient can be written 

. cosh(fi(o/2kT ) 1 f 1 

“v'VjW-^^jL^vw^ UVJ sh{ha V MllkT)n\{<o-<Ov il ,v Ji ) 2 

X + (V-Ww) + (<W + J^- <«+«W 


(40) 
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FREQUENCY u (cm' 1 ) 

Fig. 9. Comparison between the H 2 0-N 2 theoretical values and the laboratory experimental data of 
Burch et al. (1984, 1985) for T— 353 K in the spectral range 1200 cm -1 < w<2300 cm -1 . The theo- 
retical values are denoted by A and the experimental values by +, respectively. 

Because these two averaged line shape functions are common for all the bands, 
and they fall off very rapidly as their arguments increase, the continuum absorp- 
tion is mainly determined by the nearest and strongest two bands; that is the one 
on the low-frequency and the one on the high-frequency side of each window. We 
expect that the basic features of the absorption in all the windows, including the 
temperature dependence, should be, more or less, similar to each other. However, 
their magnitudes, which are roughly proportional to the total intensities of the 
nearest bands, are expected to be different, generally decreasing as the frequency 
(o increases. 


3. Computational results of the water continuum absorption 

Using the theory outlined above together with the accurate position and inten- 

sity data of the H z O transitions (Rothman et al., 1987), and the dipole moment 

of H 2 0,/r= 1.8546 D, the quadrupole moment Q= — 1.04 eal ofN 2 , and the same 
parameters of the isotropic interaction used by Rosenkranz (1985, 1987) (C/ 
k= 8.0 X 10 3 * 5 A 6 K, 5=48, and a= 3.1 3 A for H 2 0-H 2 0; and C/k=8.2x 10 5 A 6 K, 
5=9, and (7=2.94 A for H 2 0-N 2 ), we have calculated the positive and negative 
frequency resonance average line shape functions for H 2 0-H 2 0 and for H 2 0-N 2 
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Fig. 10. Comparison between the H 2 0-N 2 theoretical values and the laboratory experimental data of 
Burch et al. (1984, 1985) for T - 296 K in the spectral range 3000 cm -1 < w<4100 cm' 1 . The theo- 
retical values are denoted by A and the experimental values by + , respectively. 

as shown in Figs. 1 and 2 for T= 296K. Also shown are the corresponding results 
obtained from the theory of Rosenkranz (1985, 1987). These functions have the 
same general shape as the empirical x functions deduced from laboratory and 
atmospheric data (Clough et al., 1989) and used in many atmospheric codes. 
However, one cannot compare these functions directly because of different math- 
ematical expressions for the absorption coefficient [the use of the hyperbolic tanh 
instead of sinh, cf. Eq. (40) ], and the different treatment of the contributions 
from local lines (i.e., those within some cutoff frequency, typically ± 25 cm - 1 ). 

The corresponding theoretical absorption coefficients, a(co), for frequencies 
to 10 000 cm -1 are shown in Figs. 3 and 4, respectively, for self- and ^-broad- 
ening at T= 296 K. We note that these results vary over 5 orders of magnitude in 
this limited frequency range. Similar results have been obtained for other tem- 
peratures in the range 200 K< T< 700 K but are not shown. In order to assess the 
importance of the far-wing contribution to atmospheric measurements, we first 
compare our results to the most comprehensive set of laboratory data, i.e., those 
of Burch and co-workers (1971, 1974, 1975, 1982, 1984, 1985). 


4. Comparison between theory and experiment, preliminary atmospheric 
validations and conclusions 

The theoretical absorption coefficients for self- and N 2 -broadening are com- 
pared with the experimental results of Burch and co-workers in Figs. 5-10. From 
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Fig. 11. Comparison between the theoretical temperature dependence and the laboratory data of 
Hinderlinget al. (1987) at the 10P(20) C0 2 laser line frequency of 944. 195 cnr 1 . 

these and other comparisons ( Ma and Tipping, 1991, 1 992a, b ) , we can draw the 
following conclusions: (i) Using the different positive and negative resonance 
average line shape functions gives better results than using a single average line 
shape function; (ii) The temperature dependence for this limited range of T is 
well described by the theory (in fact, the agreement is generally better for higher 
temperatures where the measured absorption is greater due to the increased va- 
por pressure and where the corresponding experimental uncertainties are lower); 
( iii ) The agreement in the 300 < to < 1 1 00 cm ~ 1 region is very good for both self- 
and N 2 -broadening; (iv) The agreement in the 2400 <a >< 2700 cm -1 region is 
not as good, but the absolute magnitude of the continua are smaller than those in 
the 300<eu< 1 100 cm -1 region by approximately a factor of 1 00 (see Figs. 3 and 
4), and the corresponding experimental uncertainties are much larger; (v) The 
agreement for measurements made near the centers of vibration-rotational bands 
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TEMPERATURE (K) 

Fig. 1 2. Comparison between the theoretical temperature dependence and various laboratory data: O 
Loperet al. (1983); • Montgomery (1978); A Burch et al. (1971); ■ Burch data, cited by Roberts 
et al. ( 1976), all for cos 1203 cm' 1 ; O Varanasi ( 1988b) for cos 1000 cm' 1 . 

(1200-2200 cm' 1 and 3000-4200 cm’ 1 ) show greater scatter than measure- 
ments made between the bands; (vi) Most of the theoretical values are below the 
corresponding experimental measurements. We will return to this last point later. 

In order to investigate the temperature dependence in more detail over a wider 
range of T, we compare in Fig. 1 1 the theoretical calculations with the experimen- 
tal results of Hinderling et al. (1987), obtained using the 10P(20) laser line of 
C0 2 having a frequency of 944.195 cm -1 . The results are similar to those ob- 
tained using the data of Burch and co-workers; viz., the theoretical values are 
systematically lower than the experimental values, and the agreement is generally 
better at higher temperatures where the experimental uncertainties are smaller. 
Similar results obtain for comparisons made with other laser lines. In Fig. 12, we 
show a comparison for the temperature dependence of the self absorption coeffi- 
cient at frequencies near 1203 cm' 1 with experimental data obtained by a num- 
ber of researchers using different experimental techniques (Burch et al., 1971; 
Montgomery, 1978; Loperet al., 1983). Also shown are some data from Varanasi 
(1988b) for a frequency near 1000 cm -1 . It is clear that the theoretical results 
reproduce the general trend very well, but underestimate the magnitude by 30% 
or so. We note that in contrast to the empirical model of Roberts et al. (1976) or 
predictions made assuming the dimer model, the absorption coefficient shows a 
positive slope for temperatures above 400 K that is corroborated by the present 
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WAVENUMBER 

Fig. 1 3. Comparison between the Nimbus-4 IRIS satellite radiance spectrum in the 400 cm - 1 < u>< 1 600 
cm -1 spectral region and the line-by-line synthetic spectrum calculated without the inclusion of the 
far-wing water absorption; the difference is plotted in the lower part of the figure together with the 
IRIS standard deviation. 

theoretical calculations. This change of slope with temperature depends sensi- 
tively on the frequency considered. 

Bolstered by the good agreement between theory and laboratory data, we now 
discuss some preliminary validations of the theoretically calculated continua us- 
ing atmospheric measurements. In Fig. 1 3 we show the radiance of the Earth as 
measured by the Nimbus-4 IRIS spectrometer in the 400 < oj < 1 600 cm ~ 1 region. 
Also shown is the predicted synthetic spectrum obtained using a line-by-line cal- 
culation without including the effects of the water continuum (Carlson and La- 
ds, 1 992 ). For clarity, the difference between these is plotted in the lower part of 
the figure along with the IRIS standard deviation. As apparent from the figure, 
large differences are observed. In Fig. 14, we show a similar comparison except 
that the synthetic spectrum was calculated including the theoretical H 2 0-H 2 0 
continuum. The differences are dramatically reduced, especially in the important 
700-1200 cm -1 region. Many of the remaining sharp differences arise from line 
mixing (e.g. in C0 2 Q branches) or from other trace molecular species not in- 
cluded in the synthesis. We note that similar comparisons using other continua, 
for instance that of Roberts etal. (1976) or the CKD model (Clough etal., 1989) 
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Fig. 14. The same as Fig. 13 except that the synthetic spectrum was calculated by including the far- 
wing water absorption. 

would show similar reductions in the difference spectrum. This is obvious since 
the above mentioned continua are based primarily on the data of Burch et al. and 
the present theoretical values are in substantial agreement with this data in the 
spectral region covered by the IRIS instrument. Such agreement is not assured in 
other spectral regions, or more particularly, for other temperatures (Hartmann 
et al., 1993). 

As a result of the comparisons discussed above and others not shown, we can 
draw several general conclusions concerning the contributions of far-wing ab- 
sorption to the observed water continuum. The first and most important conclu- 
sion is that within the framework of the present theory and without the introduc- 
tion of arbitrary adjustable parameters, we obtain generally good agreement with 
both laboratory and atmospheric data over a wide range of frequencies and tem- 
peratures. Thus we feel that this mechanism is the major source of the observed 
continuum absorption. As noted above, while the agreement between theory and 
measurements made near the band centers is generally good, there appears to be 
a systematic underestimate of the absorption. There are two possible theoretical 
interpretations for this discrepancy: the near wings of the strong vibration-rota- 
tional lines calculated using the quasistatic approximation may differ signifi- 
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cantly from the true line shapes, or there may be another mechanism producing 
absorption near the line centers. Both of these effects are currently under inves- 
tigation. Corrections to the near-wing line shape resulting from the non-commu- 
tative term neglected in Eq. (3) and from the effects of molecular motion ne- 
glected in the quasistatic approximation would increase the near-wing absorption 
while leaving the far wings unchanged. We are also investigating the role and 
magnitude of interference effects associated with the allowed and the collision- 
induced dipoles. Such effects are known to play an important role in the spectra 
of HD (Herman et al., 1979). 

In conclusion, we note that much progress on the theoretical understanding of 
the water continuum has been made. The far-wing absorption is the primary 
mechanism responsible for the water continuum and this is described reasonably 
well within the framework of the present theory. However, improvements to the 
theory and the possible inclusion of other mechanisms to calculate the absorption 
due to the near wings of spectral lines are needed. Finally, we would like to stress 
the need for more accurate experimental data (laboratory and atmospheric). 
When such data are available, and after the refinements discussed above are made, 
one can vary the input parameters for the isotropic interaction (the Lennard- 
Jones parameters) or even implement a different isotropic interaction and see if 
the discrepancies between theory and experiment can be reduced. No such op- 
timization has been done in the present work because the differences between 
theory and experiment are well within the combined limits of errors, given the 
difficulties of the measurements and the mathematical approximations made in 
calculating the theoretical results. 
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